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Abstract Move-to-front rule is a heuristic updating a list of n items according to re- 
quests. Items are required with unknown probabilities (or popularities). The induced 
Markov chain is known to be ergodic [J. One main problem is the study of the dis- 
tribution of the search cost defined as the position of the required item. Here we first 
establish the link between two recent papers [SHI] that both extend results proved by 
Kingman [7] on the expected stationary search cost. Combining results contained in 
these papers, we obtain the limiting behavior for any moments of the stationary seach 
cost as n tends to infinity. 
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1 Introduction 

The heaps problem was first considered, in independent works, by Tsetlin [11] and 
McGabe [5]. Its basic description can be given as follows. Gonsider a collection of n 
items stored into a list or heap and each of them is identified by a label. Hence, the 
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objects can be described by the set I = {1, . . . ,n}. The probabihty that the i-th item 
is requested by a user is denoted by Pi, for i = 1, . . . ,n. Hence Pi ^ 0, for any i, and 
Pi ~ 1- each unit of time, an item is requested and it is searched for through 
the heap, starting at the top. Once it is found, it is moved to the top of the heap. 
The search cost is the position of the requested item in the heap or, equivalently, the 
number of items to be removed from the heap in order to find the requested one. In 
this setting, it might be of interest to determine the distribution of the search cost 
when the underlying Markov chain is at equilibrium. 

Kingman [7] first studied the case of random request probabilities, or random popular- 
ities. His paper develops two important cases where request probabilities are defined 
in terms of: (a) the normalized increments of a 7-stable subordinator; (b) the Dirichlet 
distribution on the simplex. The results contained therein provide an exact analytic 
evaluation of the expected search cost either for any finite n or in the limit, as the 
number of items n tends to infinity. In particular, in the case of normalized 7-stable 
request probabilities, it is found that the limiting expected search cost is finite if and 
only if 7 < 1/2. 

These results have been recently extended in two independent papers. In [8] Lijoi and 
Priinster studied the case of request probabilities derived from a normalized random 
measures with independent increments, which generalizes the result obtained by King- 
man [7] . In [5] Barrera and Paroissin studied the case of request probabilities based on 
exchangeable random partitions. 

It is to be emphasized that all previous contributions on the subject is confined to the 
determination of the first moment of the stationary search cost. Here we wish to extend 
earlier work and determine the expression of the limiting moments of any order in the 
7-stable case. In particular, it will be shown that the fc-th moment exists if and only if 
7 < + 1) which reduces to the condition provided by [7] when k — 1. See also [8]. 
The outline of the paper is as follows. In Section 2 we provide a concise introduction 
to some basic tools and notions that will be relevant for achieving the main result in 
Section 3. 



2 The 7-stable model 

Before stating and proving our result, it might be worth recalling the main ingredients 
that define the model we are going to use. As mentioned in the previous section, the 
request probabilities p,, for i = 1, . . . ,n, are going to be random. Indeed, if (iui).i>i 
is a sequence of positive independent random variables and Wn — Yl7=l '^^^ '^^^ 
define 

Pi = TTT I = 1, . . . ,n 

Wn 

Hence, (pi, . . . ,pn) is an exchangeable random partition of the unit interval. A possible 
choice is Wi := ^t- — (,ti_i where = fo < < • ■ ■ < tn = 1 and ^ = {^t '■ t G 
[0, 1]} is a subordinator that is a process with almost surely increasing paths and with 
independents increments. In this case, one can express the the Laplace transform of 
in terms of the Levy intensity v of ^. In other words 

<^,(s) := E [e-^'"-] =exp|-(t,-t,_i)^°° [l-^e-^^] v{dy)^ (1) 
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with V such that JJ^ minjl, j/} !/(dj/) < ca. According to the terminology set forth in 
[10| . (pi, . . . ,pn) defines a normahzed random measure with independent increments 
(NRMI). 

Lijoi and Priinster in [5] considered this general construction to determine an expression 
of the expected value of the search cost Sn- In the special case where 

Kdy)- r{\--i) y"^"" 7e(o,i) (2) 

they recovered an expression of the limiting expected search cost, as n tends to infinity, 
thus recovering a result proved by [?]• Note that if is as in ([2]), then (^i(s) = exp{ — (t^ — 
ti-\)s^^ for any s > 0. 

Barrera and Paroissin [3] have been able to determine an integral representation for 
the Laplace transform of the search cost Sn in terms of the Laplace transforms 
(\>i of the single random weights Wj. In doing so they rely on results proved by Fill and 
Hoist [5]. The expression they obtain is, then, used to derive a formula for the first 
two moments. From these formulas, they get an asymptotic equivalent for the Laplace 
transform of Sn and the limit of the two first moments. Only this last point needs 
the assumption that the expectation of Sn is finite. Two examples are studied: the 
case of deterministic weight and the case of gamma weight, which corresponds to the 
Dirichlet partition. Notice that, for this case, some limiting results were proved with 
an alternative way in . The limiting distribution has been also derived in the general 
iid case provided that the expectation /^^ of Wi is finite [2]. 

In the following section we will undertake the approach developed in [3] and determine 
the fc-th moment of Sn by working directly on . 



3 Moments of the stationary search cost 

The main tool we are relying on for the evaluation of EfS^] is the Laplace transform 
of Sn as displayed in theorem 2.2 of [3] and recalled here below. 

Theorem 1 For a sequence (wj)i>i of independent random variables 

n- 1^00 I POO _ _ \ 

'^S„ («) = E / / (^) n '^M,. (0 dr dt , (3) 

for all s ^ 0, where for all j G {1, . . . , n}, 

ht.s,j{'') = <t>o{r) + e-'{(j>j{r-t)^(t,j{r)) , t > 0, r > . 

Using Q, we are able to compute moments of any order of the search cost Sn- 
Before doing so we need to introduce the quantity 

rOO /"OO ^ 

X n [^o{r) + e'\4>j{r-t)-4>-j{r))] (4) 
whose values, at s = 0, will determine the moments of Sn- 
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Proposition 1 // the (pi, . . . ,pn) are determined by normalizing the increments of a 
■y-stable subordinator with ti — ti^i = 1/n in ^ for each i G {1, . . . , n}, then 

lim M,,„(0) = I (T^^ if7<^ 
I oo otherwise 

where (a)j, = r{a + k)/r{a) is the k-th ascending factorial of a. 

Proof Note first that <j)i{^) = exp{ — s'^'/n} for any s > 0. Moreover 



f-OO f-OO ^ 

Mk,n{0)= <l'^{r)\{{<l^^,{r~t)~<j,,^{r)) W ^j{r-t)drdt 

f-OO f-OO 

POO /-OO _^ _ 

X Y[i<f>^^{r)-(|'i^{r + t))dtdr 

m—1 

Taking into account the form of 0^ in the 7-stable case, one has 
k k 

m=l a,e{0,l}'="i '"=1 

^ ^_^ya,|g-<l±fl2:|a,|g-4(fc-l-|a,|) 
a(e{0,l}''"^ 

where a/ = (ai, . . . , a/_i, a;_|_i, . . . , flfc) and |a;| = X^m^/am- Summing up, in the 
7-stable case one has 



M,Ao)-^ E E E M)'"7 

X / (r + t) 



e ^ >i 



27-2^-<I±i)l(2+|aH)g-4(fc-l-|«H)didr 



S E E E (-1)'-' /'^-^^^---^ 



fc-i 



f-OO 



The change of variable {x,y) = {{r + t)'^ ,r'^) yields 



M,Mo) = ^ E E E M)'"7 



X / X 

'y 



Using formulae (3.381.6) and (7.621.3) in [6], one finds out that 

Mfe,jo) = ^ E E E (-1)'"''/ y^-Hn-\\ai\ + 2)r^^^ 

;J E E E 

E EEVir(';').i^i(2,i;i + i;i 



2+ la; 



fe k-1 , , , , 

1=1 r=0 

7fc(n-l)(n-2)---(n-fc) i . ^ , i . i 2 + r 



2F1 2,1;- + 1;1 
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Since the Gauss hypergeometric function 2F1 can be rewritten as 

.F, f 2, 1; i + 1; 1 - V iMlJ^ V(-l)^- f ' + .^ ^ 

the expression of M^. „(0) can be further simplified as follows 
Mfc,„(0) = 7fc(n-l)(n-2)...(n^fc) /fc - 1 

r=0 V ^ 



_ - 2) . . . (n - fc) /fc - A f^i-iya, 



where 

g (2)K1). M 
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A simple change of variable m — I — j leads to write aj as 

y-^ (2)m+j(l)m+j (m + j 

m— ( 



'o("^ + i)!(i + ^)™+j V 

{m + j + ly.jm + j)\ 1 



771 — •''^ 7 

oo 



(i + 1)! V- (j+2V(i + i)„(ir 



and, consequently, 



(J + 

fi + lV 1 



fc— 1 / \ CX3 

r=0 \ / j=0 ' ^ 7 •'j 

X 2Fl(j + 2,j + l,j + l + -,l) — ^ 

7 V n 

If one resorts to identity (0.154.6) in [S], it follows that 

' T), — ^ 



X 2-Fi(j + 2,j + l,j + l + i,l) ^(-1) 
Finally, using formula (0.154.5) in [6] one has 



^ 'fc-l \ /2 + r 
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1 fc!(fc- 1)! ^ , , , 1 

2i^l(fc + l,A;,fc+-,l) 



n'^-i (l + i)fe_i ' ' 7 



as n — >■ cxD. If 7 < -jriY then 2-Fl(fc + 1, fc, ^ + ^, 1) = ^ /t^iLi*)/^ (1) — ^- Otherwise the 
Gauss hypergeometric function diverges (see paragraph 9.102 in |^). After some little 
algebra the results is proved. □ 

The study of the limiting behavior of Affc_„(0) is crucial for understanding the 
limiting behavior of the moments. Indeed, 
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In particular, we have: 

E{S„) = A/i,„(0) 

E{sf,) = A/i,„(0) + 3M2^„(0) + M3,„(0) 
£(5,1) = Mi,„(0) + 7Af2,„(0) + 6A/3,„(0) + M4,„(0) 
E{S^) = Mi,„(0) + 15A^2,„(0) + 25A'/3,«(0) + 10Af4,„(0) + A4,„(0) 



B(S^) = A/i,„(0)+A^2,n(0) 



In general 



where 



E{S!^) = )A.fi,„(0) + ■ ■ ■ + a(^'^^A/4,„(0) (5) 



af-) = a^^Y^) + Za['^-i) Z = 2, . . . , fc - 1 (6) 
4^=) ^ 1 

The last recursion follows from the fact that 

Mk^nis) = -fcA4,„(s) - Mk+i,„{s) 

From proposition [T] and equation ((SJ, we have the following theorem. 

Theorem 2 // the (pi,...,p„) are determined by normalizing the increments of a 
y-stable subordinator with t^ — ti_i = 1/n in ^ for each i G {1, . . . , n}, then 



v^fc H'.r (fc) ■/■ , 1 

lim 5(5^) = <! ^'=1 TF^"' ^f'y<—i 

' oo otherwise 



The previous theorem allows to calculate all the moments of the limiting search cost 
distribution in the stable case. For example the second moment is 

7(1+7) 



lim E{S?,) = <! (i-37)(i-27) if '>'< 3 
n^oo 00 otherwise 



and the third moment 



f 7(1+57) if 7- < 1 

lim E{Sf^) = <^ (l-47)(l-37)(l-27) " T ^ 1 

n^oa QQ otherwise 
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